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Abstract. We consider a one-parameter family of invertible maps of a two-dimensional 
lattice, obtained by applying round-off to planar rotations. We let the angle of rotation 
approach tt/2, and show that the limit of vanishing discretisation is described by an integrable 
piecewise-smooth Hamiltonian flow, whereby the plane foliates into families of invariant 
polygons with an increasing number of sides. The round-off perturbation introduces KAM- 
type phenomena: a positive fraction of the unperturbed curves survives, and locally this 
fraction converges to a rational number strictly less than 1. We prove this considering 
symmetric orbits, under a condition that allows us to obtain explicit values for densities. 
Under similar assumptions, we also prove that, locally, the perturbed curves are uniformly 
distributed in phase space. The dynamics is described in terms of strip maps, a variant of 
those found in outer billiards of polygons. 



1. Introduction 

The study of near-integrable Hamiltonian dynamics on a discrete phase space presents a 
unique set of problems. On the continuum, the orbits of an integrable symplectic map are 
rotations on invariant tori, which, generically, are quasi-periodic. According to KAM theory, 
a positive fraction of these tori, identified by their frequency, will survive a sufficiently small 
smooth perturbation. The complement of KAM tori consists of a hierarchical arrangement of 
island chains and thin stochastic layers. In low-dimensions, KAM tori disconnect the space 
and hence ensure stability. 

Reproducing these structures in a discrete space is problematic, due to the lack of a frame- 
work for perturbation theory. Typically, quasi-periodic orbits do not exist. Surrogate KAM 
surfaces must thus be identified, and their evolution must be tracked, as a perturbation pa- 
rameter is varied. Even in low dimensions, these orbits need not disconnect the space, so 
their relevance to stability must be re-assessed. 

There are various approaches to space discretisation. For an algebraic system it is natural 
to replace the real or complex coordinate fields by a finite field. Because this procedure erases 
all topological information, in the discrete phase space there is no near-integrable regime at 
all, and one witnesses a discontinuous transition from integrable to non-integrable behaviour. 
This transition manifests itself probabilistically via a (conjectured) abrupt change in the 
asymptotic (large field) distribution of the periods of the orbits [131 f22] . 

Round-off in computer arithmetic brings about an equally blunt discretisation of space. 
Here a small perturbation causes the dynamics to collapse onto a discrete set. In rare cases, 
round-off fluctuations act like small amplitude noise, and give rise to Gaussian transport; the 
latter, however, wipes out all small-scale dynamical features. More commonly, the noise model 
is inappropriate, but there is no general theory to rely on. (Shadowing theory is unhelpful 
here because it requires hyperbolicity [HI Section 18.1].) In all, the literature devoted to 
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the study of deterministic (as opposed to probabilistic) manifestation of near-integrability in 
computer arithmetic is minimal [9j EJ [29] . 

A different kind of space discretisation occurs in piecewise isometric systems, as the com- 
bined effect of discontinuity and rationality. In isometries involving rational rotations, the 
space gets tessellated by polygons with algebraic number coordinates. As these polygons 
move rigidly, the phase space is discrete. Outer billiards of polygons are symplectic maps 
of this type, which feature a skeletal version of divided phase space (see [21], and references 
therein). Under appropriate rationality conditions, these systems support a countable family 
of bounding invariant sets, which prevent orbits from escaping to infinity [27\ [T5j UOj. Unlike 
KAM tori, these sets are uniformly distributed on the plane — this theme will also appear in 
the present work. The existence of unbounded orbits for some irrational parameters has been 
a significant recent advance [23] (see also [8]). 

In this paper we explore discrete near-integrability in the family of invertible lattice maps 

(1) F:Z 2 ^Z 2 (x,y) ^ ([\x\ - y, x) A = 2 cos(2vri/) 

where |_ - J is the floor function — the largest integer not exceeding its argument. If we remove 
the floor function in equation ([1]), we obtain a one-parameter family of linear maps of the 
plane, (x,y) i— > (Ax — y,x), which are linearly conjugate to rotation by the angle 2ttu. The 
floor function models the effect of round-off, pushing the image point (Xx — y, x) to the nearest 
lattice point on the left. Compared to the more straightforward rounding to the nearest lattice 
point, this procedure is arithmetically nicer, and has fewer symmetries, thereby enhancing 
the effects of the round-off. In the lattice map F, the discretisation length is fixed, and the 
limit of vanishing discretisation corresponds to motions at infinity. This asymptotic regime 
is our main interest. 

The deceptively simple model ([!]) displays a rich landscape of mathematical phenomena, 
connecting discrete dynamics and arithmetic. This model originated in dynamical system 
theory [25], [T71 [TBI 021 El 1211 HE], and was subsequently studied in number theory, within 
the context of shift radix systems [U [2] . The following unsolved question [1] [26] distills the 
difficulties encountered in the analysis of this mo del0 

Conjecture. For all real A with |A| < 2, all orbits of F are periodic. 

Due to invertibility, periodicity is equivalent to boundedness. This conjecture holds trivially 
for A = 0, ±1, where the map F is of finite order. Beyond this, the boundedness of all round- 
off orbits has been proved for only eight values of A, which correspond to the rational values 
of the rotation number v for which A is a quadratic irrational: 

±l±\/5 r r 

(2) A = ±V2, ±V3. 

(The denominator of v is 5, 10, 8, and 12, respectively.) In these cases the map F admits a 
dense and uniform embedding in a two-dimensional torus, where the round-off map extends 
continuously to a piecewise isometry, which has zero entropy (and is not ergodic). The natural 
density on the lattice Z 2 is carried into the Lebesgue measure, namely the Haar measure on 
the torus. The case A = (1 — v5) /2 was established in [17], with computer assistance. Similar 
techniques were used to extend the result to the other parameter values, but only for a set 
of initial conditions having full density [16]. The conjecture for the eight parameters ([2]) 
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was settled in [2] with an analytical proof. For any other rational value of v, there is a 
similar embedding in a piecewise isometry of a higher-dimensional torus; these systems are 
still unexplored, even in the cubic case. 

Irrational values of v bring about a different dynamics, and a different theory. The sim- 
plest cases correspond to rational values of A, and, in particular, to rational numbers whose 
denominator is the power of a prime p. In this case the map F admits a dense and uniform 
embedding in the ring Z p of p-adic integers [6]. The embedded system extends continuously to 
the composition of a full shift and an isometry (which has positive entropy), and the natural 
density in Z 2 is now carried into the Haar measure on Z p . This construct was later used to 
prove a central limit theorem for the departure of the round-off orbits from the unperturbed 
ones |28j . This phenomenon injects a probabilistic element in the determination of the period 
of the lattice orbits, highlighting the nature of the difficulties that surround conjecture [TJ 
Very recently, Akiyama and Petho |3j proved that ([I]) has infinitely many periodic orbits for 
any parameter value. 

In this work we consider a new regime, namely the limit A — > of equation (TjQ), corre- 
sponding to the rotation number v -4 1/4. This is one of five limits (the other limits being 
A — > ±1,±2) where the dynamics at the limit is trivial. After scaling, we embed the lattice 
Z 2 in R 2 , and show that there is a non-smooth integrable Hamiltonian flow (not a rota- 
tion), which represents the limiting unperturbed dynamics. Thus the limit A — > is singular. 
The parameter A acts as a perturbation parameter, and a discrete version of near-integrable 
symplectic dynamics emerges on the lattice when the perturbation is switched on. 




FIGURE 1. Some periodic orbits of the map F, for small values of the parameter A. 
Left: A = 1/25; right: A = 1/79. The boxes have side A" 1 . 

More precisely, if A is small, then the orbits lie approximately on convex polygons; the 
smaller A, the closer the approximation. The number of sides of these polygons increases 
with the distance from the origin; near the origin they are squares, while at infinity they 
approach circles. 

We introduce a piecewise smooth integrable Hamiltonian function (equation (111)) ), whose 
invariant curves are polygons, representing the limit foliation of the plane for the system 
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dU). To match Hamiltonian flow and lattice map, we exploit the fact that, for small A, the 
composite map F 4 is close to the identity. After scaling, it is possible to identify the action 
of F 4 with the unit time-advance map of the flow. The two actions agree along the sides of 
the polygons, but they differ in vanishingly small regions near the vertices. This discrepancy 
provides the perturbation mechanism. 

All integrable orbits must close after one revolution around the origin. By contrast, the 
lattice orbits need not do so, leading to a non-trivial clustering of the periods around integer 
multiples of a basic rotational period — see figure El The lowest branch of the period function 
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Figure 2. The normalised period function T\(z) for points z = (x,x), and 
A = 2 -12 . The vertical lines mark the location of critical polygons, which pass 
through lattice points. 

comprises the minimal orbits, which shadow the unperturbed orbits, and close after one revo- 
lution. For these orbits, the effects of the perturbation cancel out. (The role of cancellation of 
singularities in the existence of invariant tori of non-smooth systems was noted long ago |12j.) 
These are the orbits of the integrable system that survive the perturbation. The other orbits 
reproduce the island chains structure of a near- integrable area-preserving map, in embryonic 
form near the origin, and increasing in complexity for larger amplitudes. 

The set of invariant polygons is partitioned by critical polygons, which contain Z 2 points, 
into infinitely many polygon classes, which can be characterised arithmetically in terms of 
sums of squares. The main result of this work is that for infinitely many classes, a positive 
fraction of minimal orbits survive. We prove this by considering orbits having time-reversal 
symmetry, while the restriction to infinitely many classes — as opposed to all classes — stems 
from a coprimality condition we impose in order to achieve convergence. This phenomenon 
differs from the KAM scenario in two respects. First, as the perturbation goes to zero, the 
density of symmetric minimal orbits among all orbits does not approach unity, but attains 
a smaller rational value, which depends only on the family of polygons being considered. 
Second, the minimal orbits do not approach a Cantor-type structure, but they are uniformly 
distributed within each family. Furthermore, as the number of sides of the polygons increases 
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to infinity, the density of symmetric minimal orbits tends to zero. These results appear as 
theorems A and B, stated in section [5] after a somewhat lengthy preparation. 

Our analysis is based on the study of the first-return map to a thin strip placed along 
the symmetry axis of the map (pQ). The minimal orbits are fixed points of this map, and we 
study those lying on the symmetry line. The analysis of the return map requires tracking the 
return orbits, and this is done through repeated applications of a strip map, an acceleration 
device which exploits local integrability. This is a variant of a construct introduced for outer 
billiards of polygons (see [Ml chapter 7] , and references therein) , although in our case the strip 
map has an increasing number of components, providing a dynamics of increasing complexity. 
There is a symbolic dynamics associated with the strip map; its cylinder sets in the return 
domain are congruence classes modulo a nested sequence of two-dimensional lattices. The 
key result is that, within a polygon class, this lattice structure becomes independent of A, 
provided that A small enough. We establish this fact for infinitely many classes. This fact 
gives a 'non archimedean' character to the dynamics; the rationality of the density of minimal 
orbits, and their locally uniform distribution then follow. 

The plan of this paper is the following. In section [2] we introduce the integrable Hamil- 
tonian, characterise its invariant curves via a symbolic dynamics, and connect them to the 
arithmetical problem of sums of two squares (theorem [2]). In section [3] we switch on the 
round-off perturbation, and show that all orbits recur to a small neighbourhood of the sym- 
metry axis. Accordingly, we construct a return map of this neighbourhood, and show that 
the return orbits shadow the integrable orbits (theorem [5]) . Most proofs for this section are 
postponed to section [H Matching the symbolic dynamics of integrable and perturbed orbits 
is more delicate, requiring the exclusion of certain anomalous domains, and establishing that 
the size of these domains is negligible in the limit. This is done in section [5j where we also 
state the main results of this work, theorems A and B, and briefly discuss some experiments 
on the conditions of the statement of theorem B. In section [6] we introduce the strip map, and 
establish some of its properties (propositions [8] and [9]). In the final section we show that the 
return map commutes with a group of lattice translations. This result leads to the conclusion 
of the proof of the main theorems. 

There are many issues we haven't considered. The nature of the phase portrait at infinity, 
stability, the role played by non-symmetric orbits, the distribution of periods among the 
branches of the period function. More generally, one may consider the construct of strip 
maps for the purpose of developing a Hamiltonian perturbation theory over discrete spaces. 

These questions deserve further investigation. 

Acknowledgements. We are grateful to J A G Roberts for engaging discussions. 



Figure Q] suggests that the analysis of the limit A — > requires some scaling; equation (pQ) 
suggests that the quantity to be held constant should be Ax. Accordingly, we normalise the 
metric by introducing the scaled lattice map F\, which is conjugate to F, and acts on 
points z = \{x,y) of the scaled lattice (AZ) 2 : 



2. The integrable limit 
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The discretisation length of F\ is A. Then we define the discrete vector field, which 

A 



measures the deviation of F$ from the identity: 



(3) v: (AZ)^(AZ) 2 v(z) = F£(z) - z. 

To capture the main features of v on the scaled lattice, we introduce an auxiliary vector 
field w on the plane, given by 

(4) w:R 2 ^(AZ) 2 w(x,y) = X(2[y\ +l,-(2[x\ + 1)). 
The field w is constant on every translated unit square (called a box) 

(5) B mtTl = {(x,y) £ R 2 : [x\ = m,[y\ = n}, m,n£Z 
and we denote the value of w on B m ^ n as 

(6) w m , n = A(2n + l,-(2m + l)). 

The following proposition, whose proof we defer until section [U states that if we ignore a 
set of points of zero density, then the functions v and w agree on the lattice AZ 2 . 

Proposition 1. Let r be a positive real number. We define the set 

(7) A(r,A)={^€(AZ) 2 :||z|| 00 <r}, 
(with || (it, -y) Hoc = max(|it|, \v\)), and the ratio 

#{z£A(r,\):v(z)=w(z)} 

Then we have 

lim m(r, A) = 1. 
A— >o 

The asymptotic regime that results from replacing v by w will be referred to as the inte- 
grable limit of the system ([I]), as A — > 0. The points where the two vector fields differ have 
the property that Ax or Xy is close to an integer. The perturbation of the integrable orbits 
will take place in these small domains. 

2.1. The integrable Hamiltonian. We define the real function 

(8) P:R^R P(x) = [x\ 2 + (2[x\ + l){x} 

where {x} denotes the fractional part of x. The function P is piecewise affine, and coincides 
with the function x i— > x 2 on the integers. Thus for x G Z: 



(9) P( L xj) = Lxj 2 [^/P{x)\ = [x\. 
Using this fact, we can invert P up to sign by defining 

(10) P .K>o^K> x^ 2[ ^j +1 > 

so that (P- 1 o P)( x ) = \x\. 

We define the following Hamiltonian 

(11) &> : R 2 -)■ R &(x,y) = P{x) + P(y). 

The function & is continuous and piecewise affine. It is differentiable in IR 2 \ A, where A is 
a set of orthogonal lines given by 

(12) A = {(x, y) E R 2 : (x - [x\)(y - [y\) = 0}. 



NEAR-INTEGRABLE BEHAVIOUR 



IN A FAMILY OF DISCRETISED ROTATIONS 



7 



The set A is the boundary of the boxes B mn , defined in ([5]). The associated (scaled) Hamil- 
tonian vector field, defined for all points (x, y) G R 2 \ A, is parallel to the vector field w given 
in flU): 



(13) H r:i ^'"^?^) =w(x ' 2/) ^ eR2 \ A - 

The parameter A merely rescales the time. 

For a point z£ R 2 , we write H(z) for the level set of J 2 passing through z: 

U(z) = {w G R 2 : &>{w) = &{z)}. 

Below (theorem [2]) we shall see that these sets are polygons. The value of a polygon n(z) 
is the real number &(z), and if H(z) contains a lattice point, then we speak of a critical 
polygon and a critical value. The critical polygons form a distinguished subset of the 
plane: 

r := IJ 

zez 2 

All topological information concerning the Hamiltonian & is encoded in the partition of the 
plane generated by TU A. The elements of T act as separatrices, whose vertices belong to A. 
To characterise & arithmetically, we consider the Hamiltonian 

J2(x, y) = x 2 + y 2 

which represents the unperturbed rotations (no round-off) in the limit A — > 0. Its level 
sets are circles, and the circles containing lattice points will be called critical circles. By 
construction, the functions and J? coincide over Z 2 , and hence the value of every critical 
polygon belongs to =£?(Z 2 ), the set of non-negative integers which are representable as the 
sum of two squares. We denote this set by S '. 

A classical result, due to Fermat and Euler, states that a natural number n is a sum of 
two squares if and only if any prime congruent to 3 modulo 4 which divides n occurs with an 
even exponent in the prime factorisation of n theorem 366]). We refer to $ as the set of 
critical numbers, and use the notation 

<$ = {ei : i > 0} = {0, 1, 2, 4, 5, 8, 9, 10, 13, 16, 17, . . . }. 

There is an associated family of critical intervals, defined as 

(14) let = (ej,e i+ i). 
Let us define 

g{x) = #{e G S : e < x}. 

The following result, due to Landau and Ramanujan, gives the asymptotic behaviour of £{x) 
(see, e.g., [20]) 

(15) lim S{x) = K, 

x—>oo x 

where K is the Landau- Ramanujan constant 

1 / 1\~ 1/2 

X = T2 II =°.764.... 

¥ p prime x ' 

p=3 (mod 4) 
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Furthermore, let r(n) be the number of representations of the integer n as a sum of two 
squares. To compute r(n), we first factor n as follows 



n 



2 a n^n^ c 



where p and q are primes congruent to 1 and 3 modulo 4, respectively. (Each product is equal 
to 1 if there are no prime divisors of the corresponding type.) Then we have theorem 
278] 

(16) r(n) = 4j](& + 1) II ( ^^H ■ 

We now have the following characterisation of the invariant curves of the Hamiltonian & . 

Theorem 2. The level sets H(z) of 3? are convex polygons, invariant under the dihedral 
group Di, generated by the two orientation-reversing involutions 

(17) G : (x, y) (->• (y, x) G' : (x,y) \-+ (x,-y). 

The polygon H(z) is critical if and only if &{z) € S. The number of sides ofH(z) is equal to 

U{2n(z) + 1) &(z) S 

1 A{2n(z) + 1) - r( &{z)) &>{z) G 



(18) 



where 



(19) n(z) = \f¥{z) 

and the function r is given in U6\) . For every e£^, the critical polygon with value e intersects 
one and only one critical circle, namely that with the same value. The intersection consists 
of r(e) lattice points, and the polygon lies inside the circle. 

Proof. The symmetry properties follow from the fact that the Hamiltonian is invariant 
under the interchange of its arguments, and the function P is even: 

P{- X ) = L-xJ 2 + {-x}(2L-xJ + 1) 

(-LxJ - l) 2 - (1 - {x})(2[xj +1) x^Z 
(-fxj) 2 xGZ 

= LxJ 2 + {x}(2LxJ +1) = P(x). 

The vector field (fl~3|) is piecewise-constant, and equal to w min in the box £? TO) „ (cf. equations 
© and ©). Hence a level set H(z) is the union of segments. It is easy to verify that no three 
segments can have an end-point in common (considering end-points in the first octant will 
suffice). Thus H(z) is a polygonal curve. Equally, segments cannot intersects inside boxes, 
because they are parallel there. But a non self-intersecting symmetric polygonal must be a 
polygon. 

Next we prove convexity. Due to dihedral symmetry, if H(z) is convex within the open first 
octant < y < x, then it is piecewise convex. Thus we suppose that H(z) has an edge in the 
box B mn , where < n < m. The adjacent edge in the direction of the flow must be in one 
of the boxes 

One verifies that the three determinants 

det(w mi „, w mj „_i) det(w m>nj w m+ i )n _i) det(w min 
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are negative. This means that, in each case, at the boundary between adjacent boxes, the 
integral curve turns clockwise. So IT(z) is piecewise convex. It remains to prove that convexity 
is preserved across the boundaries of the first octant, which belong to the fixed sets FixG 
(the line x = y) and FixG' (the line y = 0) of the involutions ()17p . Indeed, Yl(z) is either 
orthogonal to FixG (in which case convexity is clearly preserved), or has a vertex (m, m) on 
it; in the latter case, the relevant determinant is det(w m _i jm , w mim _i) = —8m < 0. The 
preservation of convexity across FixG" is proved similarly, and thus Yi{z) is convex. 

The statement on the criticality of &(z) follows from the fact that, on Z 2 , we have & = =2. 

Let n(z) be a non-critical polygon. Its intersections with the x-axis have abscissas ±P~ 1 (&(z)). 
Because these numbers are not integer, using (|9|) we have that there are 2[ v /0>(z)\ +1 inte gcr 
points between them, hence as many lines orthogonal to the x-axis with integer abscissa. The 
same holds for the y-axis. It follows that YL(z) intersects A in exactly 4(2n(z) + l) points, each 
line being intersected twice. Because the vector field changes across each line, the polygon 
has 4(2n(z) + 1) vertices. 

If the polygon is critical, then we have £P(z) = e 6 $ . At each of the r(e) vertices that 
belong to Z 2 , two lines in A intersect, resulting in one fewer vertex. So r(e) vertices must be 
removed from the count. 

Next we deal with intersections of critical curves. Let us consider two arbitrary critical 
curves 

&'(x,y) = e £(x,y) = e + f e,e + /E<f. 
This system of equations yields 

{x} 2 + {y} 2 -{x}-{y} = f 
which is a circle with centre at (1/2, 1/2), and radius p, where 

(20) P 2 = f + \- 

Since we must have < {x}, {y} < 1, we find p 2 < 1/2, and since / is an integer, we obtain 
{x} = {y} = f = 0. So critical polygons and circles intersect only if they have the same value, 
and their intersection consists of lattice points. Then the number of these lattice points is 
necessarily equal to r(e). 

Finally, let the lattice point (m, n) belong to the intersection of two critical curves, and let 
w' m n be the vector field of the Hamiltonian at that point. Without loss of generality, we 
assume that (m,n) lies within the first octant. If n > 0, then the vector field of & before 
and after the vertex in the direction of the flow is equal to w m _i ; „ and w m>n _i, respectively. 
One verifies that the determinants 

det(w m _i, n , w^ n ) det(wj 7l>n , w mjTV _i) 

are negative. This means that, near this vertex, the polygon lies inside the circle. If n = 0, 
the field after the vertex is w m _i _ i, and the same result holds. 

The proof is complete. Q 

From this theorem it follows that the set Y of critical polygons partitions the plane into 
concentric domains, which we call polygon classes. Each domain contains a single critical 
circle, and has no lattice points in its interior. The values of all the polygons in a class is 
a critical interval (|14p . There is a dual arrangement for critical circles. Because counting 
critical polygons is the same as counting critical circles, the number of critical polygons (or, 
equivalently, of polygon classes) contained in a circle of radius \fx is equal to £{x\ with 
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asymptotic formula (|15p . From equation (|2Q|) . one can show that the total variation Ai2(a) 
of £2 along the polygon &{z) = a satisfies the bound 

Ai2(a) < \ 

which is strict (e.g., for a = 1). 

2.2. Symbolic dynamics of polygon classes. In theorem [2] we classified the invariant 
curves of the Hamiltonian & in terms of critical numbers. We found that the set T of critical 
polygons partitions the plane into concentric annular domains — the polygon classes. In this 
section we define a symbolic dynamics on the set of classes, which specifies the common 
itinerary of all orbits in a class, taken with respect to the lattice Z 2 . 

Suppose that the polygon Il(z) is non-critical. Then all vertices of H(z) belong to A \ Z 2 , 
where A was defined in (I12p . Let £ be a vertex. Then £ has one integer and one non-integer 
coordinate, and we let u be the value of the non-integer coordinate. We say that the vertex £ 
is of type v if [MJ = v. Then we write Vj for the type of the jth vertex, where the vertices 
of 11(2:) are enumerated according to their position in the plane, starting from the symmetry 
line FixG and proceeding clockwise. 

As the type of a vertex is defined using the modulus of the non-integer coordinate, the 
sequence of vertex types Vj reflects the eight-fold symmetry of 11(2:). Hence if the fcth vertex 
lies on the x-axis, then there are 2k — 1 vertices belonging to each quarter-turn, and the vertex 
types satisfy 

(21) Vj = v 2 h-j = V(2k-l)i+j, 1 < 3 < k, < i < 3. 

Thus it suffices to consider the vertices in the first octant, and the vertex list of H(z) is the 
sequence of vertex types 

V = . . . ,v k ). 

We note that the vertex list can be decomposed into two disjoint subsequences; those entries 
belonging to a vertex with integer x-coordinate and those belonging to a vertex with integer 
y-coordinate. These subsequences are non-decreasing and non-increasing, respectively. 

From theorem O it follows that for every e G the set of polygons II(z) with £?{z) G I e 
have the same vertex list. Let k be the number of entries in the vertex list. Since the polygon 
11(2;) is non-critical, equations (fT8|) and (fT9|) give us that 4(2[v / eJ + 1) = 4(2& — 1), and hence 

k = #V = LVeJ + 1- 

Any two polygons with the same vertex list have not only the same number of edges, but 
intersect the same collection of boxes, and have the same collection of tangent vectors. The 
critical polygons which intersect the lattice Z 2 , where the vertex list is multiply defined, form 
the boundaries between classes. The symbolic dynamics of these polygons is ambiguous, but 
this item will not be required in our analysis. 

Thus the vertex list is a function on classes, hence on £ ' . For example, the polygon class 
identified with the interval I 9 = (9, 10) (see figure [3]) has vertex list 

V(9) = (2,2,0,3). 

For each class, there are two vertex types which we can calculate explicitly; the first and the 
last. If &{z) G I e , and the polygon II(z) intersects the symmetry line FixG at some point 
(x,x) G R 2 , then by the definition (jlip of the Hamiltonian & 

&{z) = 2P(x). 
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Thus inverting P and using ([9]), it is straightforward to show that the first vertex type is 
given by 

(22) Vl = [\x\\ = [P-\&{z)/2)\ = lvQ2\. 

Similarly the last vertex type, corresponding to the vertex on the x-axis is given by 

(23) Vk =[p-\^ {z ))\ = [ y f e \. 
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Figure 3. A polygon with in the interval (9, 10) and its vertices in the 

first octant. 



3. Recurrence and return map 

(The proofs of all statements in this section are deferred until section [U) 
The lattice map F has fewer symmetries than the Hamiltonian but it is easy to verify 
that F is still reversible, being conjugate to its inverse via the involution G given in equation 

F~ 1 = GoFoG G 2 =Id. 

The scaled map F\ has the same property, and all orbits of F\ return repeatedly to a neigh- 
bourhood of the symmetry line Fix G. In this section we consider the return map associated 
to this recurrence, and identify some asymptotic properties of the first-return orbits. 
From equation ([1]), the rotation number v has the asymptotic form 

v = — arccos I - J = \ - ^- + 0(A 3 ) A -> 0. 

2vr V 2 / 4 47r 

The integer t = 4 is the zeroth-order recurrence time of orbits under F\, that is, the 
number of iterations needed for a point to return to an 0(A)-neighbourhood of its starting 
point. It turns out (see proof of proposition [T]) that the field v(z) (equation ([3])) is non-zero 
for all non-zero points z, so no orbit has period four. Accordingly, for small A > 0, we define 
the first-order recurrence time t* of the rotation to be the next time of closest approach: 

(24) £*(A) =inf{A;e N : < dtf(4i/,IN)} = - + O(l) A -> 0, 
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where djj is the Hausdorff distance, and the expression djj(x,A), with x G R, is to be 
understood as the Hausdorff distance between the sets {x} and A. 

The integer t* provides a natural recurrence timescale for F\. Let T{z) be the minimal 
period of the orbit through the point z G Z 2 . (In accordance with the periodicity conjecture, 
we assume that this function is well-defined.) Since, as A — > 0, the recurrence time t* diverges, 
the periods of the orbits will cluster around integer multiples of t*, giving rise to branches 
of the period function (figure [2]). The lowest branch corresponds to orbits that perform a 
single revolution around the origin, and their period is approximately equal to t*. The period 
function T has a normalised counterpart, given by (cf. (|24p) 

T A :(AZ) 2 ^-N T x (z) = -T(z/X). 

The values of T\ oscillate about the integers. 

We construct a Poincare return map $ on a neighbourhood of the positive half of the 
symmetry line Fix G. Let d(z) be the perpendicular distance between a point z and Fix G: 

d{z) = dn(z, Fix G). 

We define the domain X of the return map $ to be the set of points z G (AN) 2 which are 
closer to FixG than their pre- images under i 7 ^, and at least as close as their images: 

(25) X = {z G (AN) 2 : d(z) < d(F^(z)), d(z) < d(F~\z))}. 

Asymptotically, for every non-negative integer m, the set X has non-empty intersection with 
the boxes B m m: _B m _}_i m , B m m +i. The main component is in B m m: a thin strip of width 
|| w m,m|| lying parallel to the symmetry line FixG. To ensure that this component is non- 
empty, we require A < A m , where the critical parameter A m is given by 

Xm = 4(2m + 1) ' 
(See figure [4] for an illustration of this situation.) 







/ 
















''im+l 





(m,m) 

Figure 4. The main component of X in B m ^ m . The grey areas belong to the 
set A introduced in equation ([2"7]) . 
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The transit time r to the set X is well-defined for all z G (AZ) 2 apart from the fixed point 
at the origin, where we set the value of t explicitly: 



(26) r : (AZ) 2 -> N r(z) = 

Thus the first return map <3? is the function 



z = (0,0), 

mm{k G IN : Ft(z) G X} otherwise. 



We refer to the orbit of z S A up to the return time r(z) as the return orbit of z: 

r (z) = {F A fc (z) : < k < r(z)}, z G X 
We let t_ be the transit time to X under : 



t_:(AZ) 2 ^N r_(ar) 



zG{(0,0)}UX 
■ min{£; G N : F^ k (z) G X} otherwise, 

so that the return orbit for a general z G (AZ) 2 is given by: 

O r (z) = {F*(z) : r_(z) < fc < r(z)}, z G (AZ) 2 . 

To associate a unique return orbit to an integrable orbit, we define the rescaled round-off 
function R\, which rounds points on the plane down to the next lattice point: 

R\ : R 2 — > (AZ) 2 R x (w) = \R(w/\), 

where R is the integer round-off function 

i?:R 2 ^Z 2 R(u,v) = {[u\,[v\). 

For every point w G IR 2 and every 5 > 0, the set of points 

{z G (AZ) 2 : z = R x (w), < A < 6} 

that represent w on the lattice as A — > is countably infinite and unbounded. 

In the rest of this section, we state some metric properties of the return orbits, and then 
show that the return orbits shadow the integrable orbits. The proofs will be found in section 

ED 

According to proposition Q] on page [H the points of the scaled lattice (AZ) 2 at which 
integrable and non-integrable vector fields have different values are rare, as a proportion of 
lattice points. The following result shows that these points are also rare within each return 
orbit. 

Proposition 3. For any w £ R 2 , if we define the ratio 

#{z G Q T (R\(w)) : v(z)=w(z)} 

H2{W,X) = , 

#0 T (R\(w)) 

then we have 

lira H2(w, A) = 1. 

A— »0 
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To establish propositions [Q and [31 we seek to isolate the lattice points z £ (AZ) 2 where the 
discrete vector field v(z) deviates from the integrable vector field w(z). We say that a point 
z 6 (AZ) 2 is a transition point if z and its image under F$ do not belong to the same box, 
namely if 

R{Fi{z)) + R(z). 
Let A be the set of transition points. Then 



(27) 



where 



A 



U A * 



m.raSZ 



A, ; 



m,r. 




Figure 5. Structure of phase space. The boxes -Bm,™, bounded by the set A, 
include regular domains (white), where the motion is integrable. The darker 
regions comprise the set A of transition points, given in (|27p . The darkest 
domains are the components S m n of the set S C A, defined in (j36| . 



For small A, the set of transition points A consists of thin strips of lattice points arranged 
along the lines A (see figure [5]). The following key lemma states that, for sufficiently small A, 
all points z ^ (0, 0) where v(z) 7^ w(z) are transition points. 

Lemma 4. Let A(r, A) be as in equation ([?[). Then for all r > there exists A* > such 
that, for all X < X* and z E A(r, A), we have 

v(z)^w(z) => z G AU {(0,0)}. 

We conclude this section by formulating a shadowing theorem, which states that for time 
scales corresponding to a first return to the domain X, every integrable orbit has a scaled 
return orbit that shadows it. Furthermore, this scaled orbit of the round-off map converges 
to the integrable orbit in the Hausdorff metric as A — > 0. 
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Theorem 5. For any w G R 2 , let U(w) be the orbit of , and let O r {R\(w)) be the return 
orbit at the rounded lattice point. Then 

lim d H (U(w),O t (R\(w))) = 0, 

A— >0 

where du is the Hausdorff distance on R 2 . 

This result justifies the term 'integrable limit' assigned to the flow 0* . 

4. Proofs for section 02 

Proof of lemma [H Let r > be given, and let A(r, A) be as in equation ((7J). We show- 
that if A and ||z||oo are sufficiently small (and then 

v(z)^w(z) => i?(F A 4 (z)) /i?(z). 

Since 2 G A(r, A), we have 2; G -B m ,n for some |m|, |n| < [r] , where [•] is the ceiling function, 
defined by the identity \x~\ = — [— x\. Through repeated applications of F\, we have 

F x (z) =\(- y + m,x) R(F x (z)) =(-( a + l),m), 

(?R) F ^ Z) =A(-x-a-l,-y + m) R(F 2 (z)) = (-(6 + 1), -(a + 1)), 

1 > F*(z) =\(y- m -b-l,-x-a-l) R(F 3 x (z)) =( c ,-(6+l)), 

F*{z) = \(x + a + c+l,y -m-b-l) R{F^(z)) =(d,c), 

where the integers a, b, c, d are given by 

o+l = \X(y - m)], 

6+1 =rA(x + o + l)l, 

1 ' c = [X(y-m-b-l)\ 

d = [X(x + a + c+ 1)J. 

The integers a, b, c and d label the boxes in which each iterate occurs (figure [6]), and also 
give an explicit expression for the round-off term [AxJ at each step. Thus reading from the 
last of these equations, the discrete vector field v is given by 

(30) v(z) = F A 4 (z) -z = X(a + c+l,-{m + b+ 1)), 

and z is a transition point whenever at least one of the equalities d = m and c = n on the 
box labels fails. 

If the integers m, a, b, c are sufficiently small relative to the number of lattice points per 
unit length, i.e., if 

(31) max(|m|, \a + 1|, \m + b + 1\, \a + c + 1|) < 1/A, 

then the map F^ moves the point z at most one box in each of the x and y directions, so that 
the labels a, b, c and d satisfy 

(32) b,d G {m — l,m,m + 1}, a, c G {ra — l,n, n + 1}. 

Similarly, (|3ip dictates that the discrepancy between each of the pairs (b,d), (a,c) cannot be 
too large: 

(33) \b-d\,\a-c\ < 1. 
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Figure 6. The boxes in which successive iterates of z under F\ occur, for the 
typical case where m = b = d, n = a = c. 



Letting A* = 1/(2 [r] + 3), we obtain 

max(|m|, |a + l|,|m + 6+ l|,|a + c+ l|) < max(|m| + |6+l|,|a + l| + |c|) 

< max(2|m| + \b — m\,2\n\ + \a — n\ + \c— n\) + 1 

< 2\r] + 3 < 1/A*, 

so that (|32p and (|33p hold for all A < A*. Then the expression (|30p for v, combined with the 
inequality (p3|) . gives that v(z) = w(z) if and only if 

(34) m = b n = a = c. 

Suppose that z is not a transition point, so that c = n and d = m, but that v(z) / w(z), 
so that at least one of the equalities (|34|) fails. 

If a ^ n, straightforward manipulation of inequalities shows that the only combination of 
values which satisfies (f32|) is 

a = n — 1, m = 0, 6 = —1, Ay = n. 

Similarly, using also the inequality (f33|) . b ^ m implies that 

b = m — 1, c = 0, a = —1, Ax = m. 

Thus if v(z) 7^ w(z), z is a transition point unless m = a + l = fe+l = c = 0, which 
corresponds to the unique point z = (0,0). D 

We now use this result to prove propositions [Hand El given in sections [2] and [3l respectively. 

Proof of proposition [TJ From equation 0, we have that the number of lattice points 
in the set A(r, A) has the asymptotic form 

#A(r,\) = o(^) A^O, r^oo. 
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By lemma [U every non-zero point z £ A(r, A) satisfying v(z) ^ w(z) is a transition point, 
so has z € A mj „ for some m, n G Z with |m|, |n| < \r~\ + 1. Furthermore, every set A min is 
composed of two strips, each of unit length, and width approximately equal to A(2m + 1) and 
A(2n + 1), respectively. We can bound the number of lattice points in the set A m n explicitly 
by 

j2m + lj + |2n + 1| +c 

for some constant c, independent of m and n. (Indeed c = 3 is sufficient.) It follows that for 
fixed r > 0, as A — )• we have the estimate 

#{zG A(r,A) :v(*)^w(z)} 



Mi(r,A) = 1 
> 1 



#A(r, A) 
#(A(r,A)nA) 
#A(r, A) 



|m|,|n|<H+l 

= 1-0 (A 2 ) V |2m + l| + |2n + l| +c 

|m|,|n|<fr]+l 

= 1-0(A). 

Since /^i(r, A) < 1, the proof is complete. D 



Proof of proposition [3j Let w e R 2 be given, and let z = R\(w), so that C T (z) is the 
return orbit which shadows the integrable orbit Ii{w). By lemma [U for sufficiently small A, 
every non-zero point £ E T (z) with v(£) 7^ w(^) is a transition point, and the return orbit 
contains at most one transition point for every box it intersects. The number of such boxes 
is approximately equal to the number of sides (|18p of the polygon H(w). Thus the number of 
transition points in the return orbit O t (z) varies like 

#(o T (z) n A) = o (T^R) • 

as A — > and |w| — > 00. 

Now we consider the total number of points in the return orbit O t {z). As v(£) = w(£) 
whenever ^ is not a transition point, the set of points £ G O t {z) which lie in some box B m ^ n 
which does not intersect the domain X are collinear and separated by a distance |w m>n |. 
Hence their number is bounded above by 

/2 

(35) #(0^)nB m , n )<A + i» B mtn nx = (/). 

| w m,n| 

This bound requires a small alteration for the boxes B m ^ n which intersect X. The set of 
points £ G O t {z) which lie in such a box may be collinear, or may split into two collinear 
subsets, corresponding to the beginning and end of the return orbit: 

O t (z) n B m>n = {F- 4k (z) : r-(z) -i<k< r_(z)} U {Ff(z) : r{z) -j<k< t(z)}, 

for some i,j £ Z>q. Hence the cardinality of the set O t {z) n B mn is at most twice the bound 
given in (j35j) . 
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We have already noted that the number of boxes which intersect the return orbit O r [z) 



behaves like O ^(w)j. Thus the total number of points in the orbit has the following 
behaviour 

#O r (z)= #(O r (z)nB m>n ) 

< T »(tA+^ 

m,n£Z 
B m ,nnO r (i)^» 



rn.n 



o 



( y^H 



as A — > and \w\ — > oo. 

Thus for fixed io £ R 2 , we have the estimate 



H 2 {z, A) = 1 
> 1 
= 1 



#{w € Q T {z) : v(w) f w(w)} 

#(0 T (z) n A) 
#O r {z) ' 
0(1) 



0(1/ xy 

= 1-0(A). 

Since H2(r, A) < 1, the proof is complete. D 



To prove theorem [5j we bound the variation in the Hamiltonian function & along scaled 
orbits O t (z). As & is invariant along every edge of n(z), for every point z which is not a 
transition point we have 

R{z) = R{z + v{z)) => v(z) = w(z), &{z) = 0>(z + w). 

Thus & is invariant under at all points for which v(z) = w(z) and, by proposition O 
these points dominate in the integrable limit. 

PROOF of theorem El Let w e K 2 be given, and let z = R\(w), so that O t (z) is the 
return orbit which shadows the integrable orbit IL(w). For a point £ £ O t (z), we consider the 
difference between the value of ^ at £ and its value on the polygon n(io): 

As we have just seen, 2? is invariant under at any point which is not a transition point. 
For a transition point £ £ A m n , we bound the change in & under F* as A — > using the fact 
that ^ is piecewise-amne: 

PiFt (0) - ^(0 = ^(e + v(0) - ^(0 

= + O(A)) - 
= 0(A). 
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The number of transition points £ G A in the return orbit is finite and independent of A 
(approximately one for every vertex of Tl(w)), so we have 

A0 f (Ff t (z)) = O(\) 

for all points in the first quarter turn, i.e., k < r(z)/4. 

The change in & under F\ also scales like A, as if £ = X(x,y), then 

&(Fx(0) ~ = P(Hl^\ - y)) - P(Xy) 
= P(X(y-[Xx\))-P(Xy) 
= 0(X), 

where P is the piecewise-affine function defined in equation ([8]). 

It follows that A^(£) scales like A for all points in the return orbit, and hence so does the 
distance between a point and the polygon H(z): 

d H (a,n(w)) = o(x). 

Neighbouring points £,£ + v(£) in the return orbit O r (z) are also 0(A)-close as A — > 0, so 
the result follows. Q 



5. Regular domains and statement of the main theorems 

In this section we state theorems A and B, after the necessary preparation. Then we briefly 
discuss the conditions of theorem B. 

In section [3] we reduced the round-off map Fx to a return map of a thin domain X, placed 
along the symmetry axis of F\. The singular set T partitions X into sub-domains, and every 
point in each sub-domain has the same vertex list. By theorem O the return orbit O r (z) of 
a point z £ X shadows the orbit 11(2;) of the integrable Hamiltonian To ensure that the 
perturbed orbits also have the same vertex list, it is necessary to replace the elements of the 
partition of X by smaller regular domains, and then prove that, in the limit, these domains 
have full density in X. 

The positive real line is partitioned into critical intervals I e , e G S\ accordingly, we partition 
the return domain X using critical polygons. Even though, for small A, the integrable and 
non-integrable return orbits U(z) and O t (z) are close (theorem [5]) , the quantity 8? is not 
constant along O t {z) and we must ensure that perturbed orbits are confined to a single 
polygon class. This requirement gives rise to the following sequence of sets 

X e = {z£l : \/w G O r (z), &>(w) G I e } e G S. 

It remains to match the vertex list associated with I e . We start by defining the edges of 
O t (z), as the non-empty sets of the form 

Sm,n H O r (z)i m,n£Z. 

For sufficiently small A, consecutive edges of O t {z) must lie in adjacent boxes, and transitions 
between edges occur when the orbit meets the set A given in equation (|27p . Thus we call the 
set O t (z) n A the set of vertices of O t {z). By analogy with the vertices of the polygons, 
we say that the return orbit O t {z) has a vertex on x = m of type v if there exists a point 
w G O t (z) such that 

w G B mtV n F~ 4 (B m ^ ltV ) or w G An-l.v n F~ A {B m , v ). 
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Similarly for a vertex on y = n of type v. A perturbed orbit is critical if it has a vertex 
whose type is undefined, i.e., if there exists w £ O t (z) such that 

for some m,n G Z (see figure |7J). 



Figure 7. A critical vertex w £ B m>n n F A 4 (£ m ±i, n ±i). 

By excluding points whose perturbed orbit is critical, we'll construct a subset X e of X e 
with the following properties: i) for all z G A e the orbit O t (z) has the same sequence of 
vertex types as n(z); ii) the union of all sets X e still has full density in X, as A — > 0. 

We now give the construction of X e . Let the set Sc Abe given by 

(36) £ = [J £ m , n , 

m,nSZ 

where 

{zGA: || < A(||(2m + l,2n + l)|| 00 + 2)} 

and ||( 

^)^)||oo — max(|u,|, |f|). Xhe set E mB is a small domain, adjacent to the integer point 
(m, n), see figure [5j 

If z G A e for some e 6 if, we say that z is regular if two properties hold: firstly, z itself 
is not a vertex of O t (z), i.e., z £ X \ A, and secondly the orbit O t {z) does not intersect the 
set S. Points which are not regular are called irregular. Then the set X e is defined as 

(37) X e = {z G X e : G / e (A)}, 

where / e (A) C I e is the largest interval such that all points in X e are regular. 

In principle, the interval / e (A) need not be uniquely defined, and may be empty. However, 
the following proposition ensures that I e (\) is well-defined for all sufficiently small A, and 
indeed that the irregular points have zero density in X as A — > 0. 

Proposition 6. Let X e and I e be as above. Then, for all e G <§ , we have 

i im EM = i. 

A^O |I e | 

Proof. Consider z £ X such that &(z) £ I e for some e G $ . If the orbit of z strays between 
polygon classes, i.e., if z ^ X e , then we have 

3w £ O r (z) : &(w) $ I e . 

However, in the proof of theorem [5] in section U we showed that the maximum variation in 
& along an orbit O r (z) is of order A, as A — > 0: 

Vz G X, \/w £ O t (z) : 9>(w) - &>(z) = O(A). 
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Hence we have 

,38) + 0(A) -{r + ° $, P^V, 

where e' is the successor of e in the sequence S '. In both cases, &{z) is near the boundary of 

r. 

If z G X e but 2 is irregular, then either z G A or the orbit of 2 intersects the set E. If 
2 G A, then one of its coordinates must be nearly integer: 

d H (z, A) = 0(A), 

where the set A was defined in (|12|) . However, as the domain X e lies in an 0(A)-neighbourhood 
of the symmetry line FixG, it follows that both coordinates must be nearly integer, giving 

d H (z, Z 2 ) = 0(A). 

Again, 3?{z) lies in a 0(A)-neighbourhood of the boundary of I e . 
Similarly, if there is a point w G O r {z) n S, then, by construction, 

d H (w, Z 2 ) = 0(A). 

Hence we have again the estimate (|38p . 
Combining these observations, we have 

|/ e (A)l |I e \^ e (A)| , 0(A) 



\I e \ \I e \ \I e \ 

and the result follows. d 

We now give an explicit representation of X e . Let x>\ be the first entry in the vertex list of 
the critical number e. By the construction f)3Tf) of X e , we have X e C B V1V1 \ A. Hence, by 
lemma HI the discrete vector field in X e satisfies 

v(z) = w(z) = w Vim z £ X e . 

Consequently by the definition (|25p of the Poincare section X, if z = \(x, y) G X e then 

-(2ui + 1) < x- y < 2vi + 1. 

Hence the set X e is given by: 

(39) X e = {z = A(x, y) G (AZ) 2 : -(2«i + 1) < x - y < 2«i + 1, ^(af) G J e }. 

Now we show that the sequence of sets X e fulfil their objective, which was to exclude all 
points z £ X whose perturbed orbit is critical in the sense defined above. 

Proposition 7. If z G X e for some e G $ , then the perturbed orbit of z is not critical. 

Proof. Suppose that z G X with 8?{z) G I e is critical, i.e., there exists w G O r (z) such that 

w G B mtU and (w) = w + v(w) G S m ±i )Tl ±i 

for some m,n G Z. We will show that z ^ X e . For simplicity, we assume that m and n are 
both non-negative, so that by the orientation of the vector field in the first quadrant: 

F*(w) = w + v(w) G S m+ i )H _i. 
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Recalling the proof of lemma 01 the expression (|30p for the perturbed vector field v at the 
point w = \(x, y) implies that 

(40) v(w) = A(o + c + l,-(m + 6 + l)), 

where the integers a,b,c are given by ([29]) , By assumption, w + v(w) E 2? m+ i in _i, which 
implies that 

c = n — 1 and d = m + 1. 
It follows that the difference between the values Ax, Ay and the integers m + 1, n, respectively, 
is bounded according to 

-A(a + n) = -A(a + c + 1) < Xx - (m + 1) < 

< Ay - n < A(m + 6+1). 
Combining this observation with the bounds (1321) on a and b gives 
||u> — (m + 1, n)||oo <A max(|a + n\, \ m + 6 + 1|) 

< Amax(|2n + 1| + \n — a\ + 1, \2m + 1| + \m — b\) 

< A(||2m + l^n + lIU + 2). 

Hence w £ £ and z is irregular, so z ^ X e . The cases where m or n are negative proceed 
similarly. Q 

Below, in equation (f56|) of section [71 we shall define a sequence of lattices L e C (AZ) 2 , 
e G (f, independent of A up to scaling, such that within the domain X e , the return map <!> is 
equivariant under the group of translations generated by L e . Formally, 

Theorem A. For every e E <§ , and all sufficiently small A, the map $ commutes with 
translations by the elements of L e on the domain X e 

(41) VlGL c , Vz,z + IeX e , $(z + 1) = &(z) + I (modw^J. 

There is a critical value of A, depending on e, above which the statement of the theorem 
is empty, as X e is insufficiently populated for a pair of points z, z + I G X e to exist. We 
write (mod a) for some vector a to denote congruence modulo the one-dimensional module 
(a) generated by a. The congruence under the local integrable vector field w Vl iVl in equation 
(jHJ is necessary for the case that $(z) + 1 ^ X. 

Furthermore, we define the fraction of symmetric, minimal orbits in X e 

#{z£X-:G{0{z)) = Q{z)Mz) = z} 

0{e,X) = = , 

#X e 

and prove the following result on the persistence of such orbits in the limit A — > 0. 

Theorem B. Let e G $ , and let (vi, . . . ,Vk) be the vertex list of the corresponding polygon 
class. If 2v\ + 1 or 2vk + 1 is coprime to 2vj + 1 for all other vertex types Vj, i.e., if 

(42) gcd(2vi + 1, 2vj + 1) = 1 i G {1, k}, 1 < j < k, Vj ^ Vi 

then, for sufficiently small A, the number of symmetric fixed points of in X e modulo L e is 
independent of A. Thus the asymptotic density of symmetric fixed points in X e converges, 
and its value is given by 

(43) lim 6(e, A) = — == . 

(2Lv^J+i)(2Lv^72J + i) 
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As in theorem A, the smallness of A serves only to ensure that X e is sufficiently populated 
for all congruence classes modulo L e to be represented. 

The condition (|42|) on the orbit code is clearly satisfied for infinitely many critical numbers 
e, e.g., those for which either 2 [v^J + 1 or 2|_y / e/2j + 1 is a prime number. The first violation 
occurs at e = 49. Numerical experiments show that the density of values of e for which (|42p 
holds decays very slowly, reaching 1/2 for e ~ 500, 000. 

The stated condition on the orbit code is actually stronger than that we require in the 
proof. This was done to simplify the formulation of the theorem. We remark that the weaker 
condition is still not necessary for the validity of the density expression (I43p . At the same 
time, there are values of e for which the density of symmetric minimal orbits deviates from 
the given formula, and convergence is not guaranteed. Our numerical experiments show that 
these deviations are small, and don't seem connected to new dynamical phenomena. More 
significant are the fluctuations in the density of non-symmetric orbits. Its dependence on e 
is considerably less predictable than for symmetric orbits, see figure El 



Density of Fixed Points 
Density ot Fixed Points on FixG 
1/(2l«J+1) 



0.02 ■ 



' 1 ' 1 1 1 1 

20 40 60 80 100 120 

e 

Figure 8. The density of minimal orbits on the symmetry line, as a function 
of the critical number e. The solid line represents the estimate (I43|) . The 
scattered points correspond to the density of all minimal orbits, symmetric 
and non-symmetric. 



6. The strip map VP 

In section [3l we saw that all points z G (AZ) 2 where the discrete vector field v(z) deviates 
from the Hamiltonian vector field w(z) lie in the set of transition points A, defined in (I27p . 
In order to study the dynamics at these points, where the perturbations from the integrable 
limit occur, we define a transit map ^ to A which we call the strip map: 

f:(AZ) 2 ^A tf(z) =F* t(z \z), 
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where, as for the transit time r to X defined in (|26p . the transit time t to A is well-defined 
for all points excluding the fixed point at the origin: 



t(z) 



z = (0,0) 

min{i 6 N : F^(z) £ A} otherwise. 

By abuse of notation, we define to be the transit map to A under i*T . Note that 
is the inverse of \F only on A. 

If z £ B mt n \ A for some m, n 6 Z, then lemma [41 (page [T4")) implies that \F(z) satisfies 

(44) ty(z) = z + t(z)w m , n , 

where w m]n is the value of the Hamiltonian vector field w in the box B m ^ n . If z £ A m;n , then 
we may have v(z) ^ w(z), so the expression becomes 

(45) $>(z) = z + v(z) + (t(z)-l)w m , n . 

In the previous section, we identified the set O r (z) n A as the set of vertices of the perturbed 
orbit O r (z). Thus, within each quarter-turn, the strip map VP represents transit to the next 
vertex. For 1 < j < k, where k is the length of the vertex list at z, we say that the orbit 
O t {z) meets the jth vertex at the point *$> J (z) £ A. For z £ X regular, the polygon II(z) 
and the return orbit O r {z) are non-critical, by construction, and the number of sides of each 
is given by equation (fl8|) . Thus the full set of vertices of O t {z) is given by 

3 2fc-l 

o T (z)nA={J J {(*'<>**)(*)}. 

i=0 j=l 

Recall that the vertices of a polygon (or orbit) are numbered in the clockwise direction — the 
orientation of the integrable vector field w. Hence the first 2k — 1 vertices (those lying in 
the first quarter-turn) are given by (^/ J (z)), 1 < j < 2k — 1. The action of F\ moves points 
from one quadrant to the next in the opposing (anti-clockwise) direction, so that the vertices 
((\F ? o F)(z)), 1 < j < 2k — 1 are the last 2k — 1 vertices. Thus the following proposition is a 
simple consequence of the number of vertices of a given polygon class. 

Proposition 8. Let e be a critical number, and let k be the length of the vertex list of the 
corresponding polygon class. Then the return map <3? on X e is related to VP via 

(46) $(z) = (V 2k oF x )(z) (mod w^), z £ X e 

where v\ is the type of the first vertex and w Vl>Vl is the value of the integrable vector field w 
at z. 

Proof. Let z £ X e . By the preceding discussion, the last vertex in the orbit O t (z) is given 
by 



i^- 1 o F x )(z) £ A V1 



Vl ■ 



The point $>(z) satisfies $>(z) £ B Vl}V1 \ A. Using the expression (j44l) for ^ applied to $>(z), 
we have 



(^ 2k oF x ){z) ee $(z) (modw„, 



VI , 



as required. D 

For z £ X regular, we use the vertices (vF 7 (z))i<j<2k-i m the first quarter-turn to define 
a sequence of natural numbers o~(z) called the orbit code of z, which encapsulates how the 
asymptotic orbit O t (z) deviates from IT (z). 
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Suppose the jth vertex of n(z) is a vertex of type Vj lying on y = n, and the orbit O r (z) 
meets its corresponding vertex at ^ 3 {z). We define the pair (xj,yj) via 



(47) ¥(z) = A 



A 



+ x 
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where Xj > 0, and \yj\, which is (essentially) the number of lattice points between *fl 3 (z) and 
the line y = n, is small relative to 1/A. Using similar arguments to those in the proof of 
proposition [7] one can show that yj satisfies 

-(2vj + 1) < yj < or < yj < 2vj + 1, 

depending whether the integrable vector field is oriented in the positive or negative y-direction. 
In both cases, the possible values of yj form a complete set of residues modulo 2vj + 1. Hence 
the jth element aj of the orbit code a{z) is defined to be the unique residue aj G {0, 1, . . . , 2vj} 
which is congruent to yj: 

(48) aj = yj (mod 2vj + 1) . 

We call y the integer coordinate of the vertex and x the non-integer coordinate. Simi- 
larly, if the jth vertex lies on x = m, then the jth element aj of the orbit code is defined to 
be the residue congruent to Xj modulo 2vj + 1 . In this case x is the integer coordinate and y 
is the non-integer coordinate. 

For all vertices in the first quadrant, the fact that orbits progress clockwise under the action 
of means that yj will be non-negative wherever y is the integer coordinate, and Xj will be 
negative wherever x is the integer coordinate: 

(49) - (2vj + 1) < Xj < or < yj < 2 Vj + 1. 
Thus the value of aj is given explicitly by 

(50) aj = Xj + 2vj + 1 or aj = yj, 
respectively. 

In addition to the values aj for 1 < j < 2k — 1 we consider cr_i, which corresponds to the 
last vertex before the symmetry line, i.e., to the point ^>~ 1 (z). Thus the orbit code of z is a 
sequence a(z) = (a-i,a%, . . . , c^fc-i), such that 

< cr_! < 2vx + 1, 

< aj < 2vj + 1, 1 < j < 2k - 1, 

where the Vj are the vertex types. 

The regularity of z ensures that the discrete vector field v matches the Hamiltonian vector 
field w in the integer coordinate at ^ J (z). The possible discrepancy in the non-integer 
coordinate is determined by the value of aj. 

Proposition 9. Let e be a critical number and let k be the length of the vertex list of the 
corresponding polygon class. For any z £ X e and any j £ { — 1, 1,2, ... ,2k — 1}, let m, n be 
such that \E ,J (z) G A mjra . Then the discrete vector field at the jth vertex of O t {z) is given by 

v(^(z)) = w mj „ + Xej(aj)e, 

where ej is a function of the jth entry aj of the orbit code a{z) and e is the unit vector in the 
direction of the non-integer coordinate of the vertex. 
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Proof. If z G X is regular, then by proposition [7J the perturbed orbit O(z) is not critical. 
Thus for any vertex w, which, by construction, satisfies 

W £ k m>n Fx(w) = W + v(u>) G B m>n 

for some m,n G Z, we must have either w G £? m ,n±l or u> G B m ±l,n- For definiteness 
we suppose that if G B m ^ n+ i, so that the vertex to lies on y = n + 1. The cases where 
w G B m>n -i or w G B m ±i^ n are similar. 

Now the proof proceeds very much as that of proposition [7J The perturbed vector field 
v(w), with R(w) = (m, n + 1), is given by equation ppj) . with a, 6, c, cZ as in (f29]h In this case, 
R(F^(w)) = (m,n) implies that 

c = n and d = m, 
and according to (j32[) . the remaining integers a and 6 satisfy 

a G {n, n + 1}, b = m. 

Thus we have 

(51) v(w) = A(n + a + l,-(2m + l)) 

= w mj „ + A(a - n)e, 

where e = A(1,0) is the unit vector in the x-direction, the non-integer coordinate direction 
of the vertex. 

If w = tyi(z), then the coefficient of the difference between v(w) and w mj „ in the x-direction 
is given by 

€j = a — n 

= \X(y-m)] -(n + 1) 

1 Xy — (n + 1) > Am 
otherwise. 



As in equation (j47j, we write 



y 



n+l 



A 

where yj satisfies the second inequality in (jMJ). Hence, by (|50|) . we have = (Tj, and the 
function = €j((Tj) is given by 

1 <jj > m 
otherwise 

which completes the proof. Q 

Note that the function ej depends on j via m. In what follows we shall write ej, omitting 
the argument. 

Applying proposition [9] to equation (jl5j) . we have that if ^ 3 (z) G A mj „, then the transit 
between vertices satisfies 

(52) ^ +1 {z) = &(z) + Ae,e + tw m , n , 

where t = t^^z)) is the transit time. Hence we think of an orbit as moving according to the 
integrable vector field at all points except the vertices, where there is a mismatch between 
integrable and non-integrable dynamics, and points are given a small 'kick' in the non-integer 
coordinate direction. 
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7. Lattice structure and proof of main theorems 
In this section we prove theorems A and B, stated in section [5j 

For egl, suppose the vertex list V(e) = (v%, . . . , vt) contains I distinct entries. We define 
the sequence (t(j))i<j<z suc h that the t(j)th entry in the vertex list is the jth distinct entry. 
Since all repeated entries are consecutive, it follows that the vertex list has the form 

(53) V(e) = (w t( i),...,t> t(1 ),t; t (2),...,u t (2), ••• > v l(1), ■ ■ ■ > v l(1))i 
with = v\ and v L n\ = v^. We define the vector L = L(e, A) as: 

(54) L= Xq (1,1), 
where the natural number q = q(e) is defined as follows 

(55) q = lcm((2u t(1) + l) 2 , (2v l{1) + l)(2u t(2) + 1), . . . , (2« t( ,_i) + l)(2u t(0 + 1)). 

Here the least common multiple runs over (2v± + l) 2 and all products of the form (2vj + 
l)(2fj + i + 1), where Vj and Vj + \ are consecutive, distinct vertex types. Finally, the lattice 
L e = L e (A) C (AZ) 2 of theorem A is given by 

(56) L e = /l, \ (L - w Vl , v 

where (• • • ) denotes the Z-module generated by a set of vectors, and the vector w Vl>Vl given 
by © is the Hamiltonian vector field w in the domain X e . 




Figure 9. The lattice L e . 
Any z, z € X e which are congruent modulo L e are related by: 

z = z + - ((2a + 6)L - 6w„ 1)t)1 ) , 

where a, b G Z are the coordinates of z — z relative to the module basis. We note that the 
vector L is parallel to the symmetry line FixG, and hence parallel to the strip X e , whereas 



28 



HEATHER REEVE-BLACK AND FRANCO VIVALDI 



the vector w VlVl is perpendicular to it. It follows that both a and b are determined uniquely 
by the value of the coefficient 2a + b, because if z = A(x, y), then 

x>y => be {0,1}, 
x < y bE {-1,0}. 

The point z itself corresponds to a = b = 0. 

We prove theorems A & B via several lemmata. The first and most significant step is to 
show that the orbit codes o~(z) of points z € X e are in one-to-one correspondence with the 
equivalence classes of X e modulo L e . We do this by constructing a sequence of nested lattices 
whose congruence classes are the cylinder sets of the orbit code. 

We define recursively a finite integer sequence (qj), j = 1, ... ,2k — 1, as follows: 

qx = (2 Vl + l) 2 

(57) Q={ qj - 1 if ^ = ^-i ->! 
H3 l k lcm((2t; J - + l)(2u i _i + l), ?J -_i) Mv^Vj-x 

Then we let 

(58) p J =q J /(2v j + l) j = l,..., 2k - 1. 

By construction, pj is also an integer. After defining the associated sequence of vectors 

(59) L 3 . = ^L_(1,1), 
we let the lattices Lj be the Z-modules with basis 

(60) Lj = (jjj, ~ (Lj — w Vl)V 
By construction 

^ fc _ 1 cq,_ 2 c...cqc(AZ) 2 . 
We claim that for all 1 < j < 2k — 1, the closed form expression for qj is given by 

(61) Qj = lcm((2« t(1) + l) 2 , (2v L(l) + l)(2u t(2) + 1), . . . , {2v i{l _ l) + l)(2« t(f) + 1)), 

where i is the number of distinct entries in the list {v\,V2, ■ ■ ■ ,Vj). That the lowest common 
multiple (|6lT) runs over all products (2vj + l)(2vj + i + 1) of consecutive, distinct vertex types 
follows from the form (153j) of the vertex list and the symmetry (|21l) of the vertex types. 
Furthermore, since all distinct vertex types occur within the first k vertex types, the expression 
(I6ip implies that the sequence (qj) is eventually stationary: 

(62) qj = q, Lj = L e , k < j < 2k - 1, 

where q and L e are given by equations (f55|) and (|56]l . 

For given e, the following result details the role of the as cylinder sets of the orbit code. 
Applying the result for j = 2k — 1, along with the observation (|62p . implies that two points 
share the same orbit code if and only if they are congruent modulo L e . 

Lemma 10. Let e be a critical number, let k be the length of the vertex list of the corresponding 
polygon class, and let pj and Lj be as above. For any 1 < j < 2k — 1 and all z,z € X e , the 
following three statements are equivalent: 

(i) the orbit codes of z and z match up to the jth entry, 



NEAR-INTEGRABLE BEHAVIOUR IN A FAMILY OF DISCRETISED ROTATIONS 29 

(ii) z and z are congruent modulo Lj, 

(Hi) the points ^ 3 (z) and ^ 3 (z) are congruent modulo Xpje, where e is the unit vector in the 
direction of the non-integer coordinate of the jth vertex. 

Proof. Let e G £, let be as above, and let z,z£ X e with orbit codes (er_i, <ti, ■ ■ ■ , o~2k-i), 
and (Ti, . . . , a2fc-i), respectively. We proceed by induction on j, with two induction 

hypotheses. Firstly we suppose that (i) is equivalent to (ii), so that for any 1 < j < 2k — 1: 

(HI) (<T_i, o-i, . . . ,0-j) = (<7_i,cti, . . . ,5-j) 44> z^z (mod L|) . 

Thus any such z is related to z via 

(63) 5 = z + - ((2a + b)Lj - bw Vuvl ) , 

for a G Z and 6 € {0,1} or 6 € {0,-1}, as appropriate. Secondly, we suppose that (ii) is 
equivalent to (Hi). In particular, 

(H2) z = z + -((2a + 6)Lj -bw Vl}V1 ) ¥ (z) = ¥ (z) + X(2a + b)pje, 

where e is the unit vector in the direction of the non-integer coordinate of that vertex. 

We begin with the base case j = 1. Suppose that the first vertex of a polygon in class e lies 
on y = vi, so that y is its integer coordinate (if x is the integer coordinate, then the analysis 
is identical). By symmetry, the previous vertex lies on x = v\ and its integer coordinate is 
x. Using the properties of ^ given in equations (@J|) and ([15]). applied to z G B VljVl \ A and 
(z) G A„ 1)Vl respectively, we have 

(64) *(z) = z (modw„ liBl ), 

(65) ^- 1 (z) + V (^- 1 (z)) = z (modw^). 
Furthermore by proposition [9[ 

v(* _1 (z)) = vr Vl}Vl + Ae_ie, 

where e = (0, 1) is the non- integer coordinate vector for the (— l)th vertex. Thus if z = 
A(|Vl/A] + x, |Vl/A] + y), by the definition ([15]) of the orbit code, the x- and y-components 
of equations ([65]) and ([6l]) . respectively, give us that the first two entries in the orbit code 
o~(z) satisfy 

x = (mod 2v\ + 1) , 

y = o\ (mod 2v\ + 1) . 

It follows that z, z G X e share the partial code (c-i, cri) if and only if 

z = z (mod (A(2ui + 1)Z) 2 ) . 
The lattice Lf is given by (cf. (160]) and <\60h ) 

L? = ^ L i, ^ ( L i - w t>i,fi)^ , 
where Li = Api(l, 1), p\ = 2ui + 1 and 

^ (Li - w VltVl ) = ~(pi -Pi,pi +Pi) = Apie. 
Thus Lf = (A(2u 1 -(- 1)Z) 2 and the first hypothesis holds. 
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Now let z, z G X e satisfy with j = 1. If $(z) = F^(z) = z + tw VljV1 , where t G IN is 
the transit time to A, then the identities 

5 + (t + a + b)w VljVl = + ^(2a + 6) (Li + w» lit)1 ) 
= + X(2a + b)pie 

with e = (1, 0) show that z has transit time t+a+b , and therefore VP (5) = ^(z) + X(2a+b)pie, 
as required (see figure [T0|) . This completes the basis for induction. 




Figure 10. The points z, z, VP(z) and *ff(z). 

Now we proceed with the inductive step. Let (|H1|) and (|H2j) hold for some j > 1. Then 
z and 5 are related as in equation ([63]) . for some a, b. We think of <x,+i, the (j + l)th entry 
of the orbit code of z, as a function of (a, 6). We suppose that the jth vertex lies on y = n 
for some n G Z (again the case in which the vertex lies on x = m is identical). Let the pair 
(xj,yj) be defined from 9 3 (z) via equation ([37]). Similarly, ^(5) defines the pair (xj,ijj). 

By ([HID, satisfies 

= + A(2a + &)pj-x, 

where x = e = (1,0). Combining this expression with equation (|52|) . applied to *$> J (z) G 
Av^n—x, we obtain 

& +1 (z) = &(z) + Ae jX + tw VJ>nr -i 
(66) = W(z) + A(2a + b) Pj x + Xejx + tw Wj , n _i. 

where t is the transit time of z to A. 
There are now two cases to consider. 

Case 1: Vj = Vj+i. 

In this case the jth. and (j + l)th vertices lie on parallel lines, which we take to be y = n 
and y = n — 1, so ^ 3+1 (z) is given by 



NEAR-INTEGRABLE BEHAVIOUR 



IN A FAMILY OF DISCRETISED ROTATIONS 



31 



and similarly for tyi +1 (z). According to the definitions ([58]) and (f60|) . we have pj = Pj+i and 
Lj = L| +1 . Thus, to show that (jHip continues to hold, we need to show that cfj+i = Cj+i for 
all (a, 6). Similarly we need to show that the vector ^>i +1 (z) — (z) is equal to the vector 
&(z) - &(z) of hypothesis (|H2|) . 

Because y is the integer coordinate of both jth and (j + l)th vertices, the transit time is 
the same for the orbits of z and z. Therefore equations (]52p and (I66p with £ = t give us that 

qji +1 (z) = W(z) + A(2a + b)p,x + Ae^-x + iw„ jin _i, 
= V j+1 (z) + \(2a + b) Pj x, 

and the second hypothesis (IH2j) remains satisfied. Furthermore, (S) and have 
the same integer (y) coordinate. It follows that, by the definition (j48]) of the orbit code, 
= Vj+i and (jHip is also satisfied. 
By the y-component of (|66p . the value of 0j+i is determined explicitly by the congruences 



n — 1 




~n 




A 




A 





(mod 2vj + 1) . 



Equation (j67|) shows that, if Vj = fj+i, then there is a map Oj i-> <Tj+i, which is a bijection 
of a set of congruence classes. 



Case 2: Vj ^ Vj+i ■ 

In this case the jth and (j + l)th vertices lie on perpendicular lines. We take these to be 
the lines y = n and x = Vj + 1, respectively, so that Vj + \ = n — 1 and is given by 



W +1 (z) = A 



A 



+ a; 



'i+i, 



n — 1 
A 



(If x is the integer co-ordinate, then the analysis is identical.) We shall demonstrate the form 
of Lj +1 by identifying those pairs (a, b) for which Oj + i = <7j+i. 

Taking the x-coordinate of equation ()66|) . and recalling the explicit form (|50p of the orbit 
code, we see that cr^+i is determined by 



(68) 



A 



+ xj + (2a + 6)pj- + 6j + t(2u i+ i + 1) 



VJ + 1 

A 



f2t> 



'i+i 



+ D 



We think of this as an integer equation of the form A(2a + b) + Bt = C , which has solutions 
2a + b S Z and i G hi for some given value of if and only if 



C 



A 



+ a j+ i 



(2v 



'i+i 



+ D 



A 



is sufficiently large and C = (mod gcd(A, £?)), i.e., if A is sufficiently small and satisfies 
the congruence 



(69) 



A 



+ Xj + Cj 



A 



(mod gcd(p 3 -,2«j + i + 1)) . 



To find the lattice L| +1 , we need to solve this equation in the case (Tj+i = crj+i- 

By assumption, the point z, given by the module coordinates a = b = 0, corresponds to 
the solution 2a + b = 0, t = t, for some transit time i G INI. Hence the general solution of ([68]) 



32 



HEATHER REEVE-BLACK AND FRANCO VIVALDI 



is given by 

(70) t = t - Pj 



gcd(pj,2v j+1 + 1) 

(71) 2a + b 



2v j+ i + 1 



gcd(pj,2v j+ i + 1)' 

for s £ Z. The second of these equations implies that s must have the same parity as 2a + 6, 
so we can write s = 2a + b, where a £ Z and 6 6 {0, ±1} for an appropriate choice of sign. 
Substituting this expression into equation ([63jh the points z for which = (Jj+i are given 
by 

1 / 2v j+1 + 1 

2 V s gcd( Pj ,2t; i+1 + l) 



V gcd( Pj ,2v i+ i + 1) 
z + - ((25 + 6)Lj + i - bw Vl)t)1 ) . 



The last equality is justified by the identities 

2^+1 + 1 L = (2^- + i + l)(2r; j + 1) Mj_ Q ^ 

gcd^-, 2v j+1 + 1) 3 gcd((2u J - + l) Pi , (2« i+1 + l)(2 Wj + 1)) 2vi + 1 1 ' J 



Mj+i 



(l,l)=L j+1 



2ua + 1 

where we have used the relationship lcm(a, b) = ab/ gcd(a, b). Therefore the first hypothesis 
(IHlj) remains satisfied. 

Substituting the general solution (i70j) and (fTTI) into equation (|66|) with x = e = (1, 0), and 
using equation (f52|) . we find 



^' +1 (5) = *i( z ) + A(2a + 6) Pj x + Aejx + tw % - i()i+1 

- + gcd(,^ +1 + l) ^ + ^ " W — ) 

= ^+ 1 (z) + \(2a + b) Pj+l y 

where y = (0, 1). Thus the points where these z meet the (J + l)th vertex share the same 
integer coordinate. So hypothesis (|H2p also remains satisfied, completing the induction. 

Thus hypotheses (|H1|) and (|H2|) hold for all 1 < j < 2k — 1, and the equivalence of the 
three statements follows. Q 



Corollary 11. Let e be a critical number, and let k be the length of the vertex list of the 
corresponding polygon class, and let j be in the range 1 < j < 2k — 1. Then two points z 
and z in X e have the same orbit code if and only if the points $ J (z) and & 3 (z) are congruent 
modulo Xq/ (2vj + 1) e, where e is the unit vector in the direction of the non-integer coordinate 
of the j th vertex. 

PROOF. Recall that for all 1 < j < 2k — 1: 



L e C Lj. 
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Thus any two points which are congruent modulo L e are also congruent modulo L|. In 
particular, if 

z = z + - ((2a + b)L - bw VuVl ) , 

= z + - ( (2a + b)—Lj - bw vltV1 
I \ qj 

then by the hypothesis (|H2|) of lemma 1101 we have: 

= + A ^(2a + 6)^ e 

= ^(z) + A(2a + 6)-^-e 

as required. D 

Lemma [TOl shows the equivalence between orbit codes and congruence classes of L e . To 
complete the proof of theorem A, we show that the orbit code cr(z) determines uniquely the 
behaviour of z under the return map $. 

Proof of theorem A 

Consider two points z, z + 1 £ X e for some I S L 6 given by 

i = -((2o + 6)L-6w Vx>Vl ). 

By corollary 111! these two points have the same orbit code and reach the (2k — l)th vertex at 
the points ty 2k ~ 1 (z), 1 $> 2k ~ 1 (z + l) G A vl r Vl+ i\, which are congruent modulo Ag/(2 , U2fe-i + l)e, 
where e is the unit vector in the non-integer direction. In particular, y i/ 2k ~ 1 (z) and ^ 2 
are related via 

m 2k -\ z + l) = m 2k -\ z ) + A(2a + 6)p2 fc _ie, 
(72) =m 2k -\z) + \{2a + b)- 



(2«i + l) 

where we have replaced U2fe-i by ui using the symmetry (|2ip of the vertex types. We will 
show that the points where they reach the last vertex are related by a similar equation: 

(73) ^ 2k - 1 o F x )(z + l) = i^- 1 o F x )(z) + A(2a + b) g e 1 , 

(2«i + 1) 

where the unit vector e^, the non- integer direction of the last vertex, is perpendicular to e. 

The last vertex of the return orbit O t (z) lies in the set A„ li „ 1 , so must be close to the image 
under F\ of the (2k — l)th vertex (see figure [Tl]) . If the (2k — l)th vertex lies on the line 
x = Vi + 1, it is a simple exercise to show that these two points are in fact equal, i.e., that 
(\j/ 2fc -! o F\)(z) = (F\ o ^ 2k ~ l )(z) for any z £ X e . We consider the less obvious case in which 
the (2k — l)th vertex lies on the line y = —v\ and the non-integer direction is e = (1,0). By 
the orientation of the vector field in the fourth quadrant, the orbit of the point z reaches this 
vertex at the point ^ 2k ~ 1 (z) given by: 



(74) * 2 *" 1 (z) = A 

where X2k-i > and < y2k-i < % v i + 1- 



vi 
A 



+ V2k-1 
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A U1 , V1 




(v,,-(I, 4-l» 



Figure 11. The (2k — l)th and last vertices, joined by the action of Fx- 



Applying F\ to (fT2j) and substituting the expression ([TI|) . we get: 



(F x o + = (Fx o ^-^(z) + A(2a + 6) 



(2^ fc + 1) p 
(2vi + 1) 













^1 - 


A 


- 2/2/5-1! 


A 



(2v k + 1) P 
' (2ui + l) 



where the non-integer direction of the last vertex is e 1 - = (0, 1). By equation (|49p . if the first 
component of this point satisfies 



A 



Vx 

A 



<0, 



-(2ui + 1) < Ui - 

then the point is the last vertex on x = v\ that we seek for all (a, b): 
(75) (tf 2 *- 1 o F A )(* + = (Fx o ^ 2k ~ l )(z + /). 



As is non- negative, U2k-i satisfies < U2k-i < 2v\ + 1, and the absolute value of the two 
ceiling functions differ by at most one, the lower bound of this inequality cannot be violated. 
For sufficiently small j/2fc-i the upper bound may be violated, in which case (Fx o ^ 2k ~ l )(z + 
I) G B VljVl and we apply F^ 4 to find: 



(tf 2 *- 1 o + I) = (F~ 3 o ^ 2fc -!)(z + /) 

-,2k-l\ 



(76) 



(F A o ^~')(z + 1)- v((F~ s o ^-^(z + I)) 
(Fxo* 2k - 1 )(z + l)-w VuVl -\ee ± , 



where the error term e is independent of (a, b) by proposition In both cases ([75]) and ([75]) 
the relationship ([75]) follows. 
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Using the property (|45p , and substituting the expression ([73]) into P6|) , we obtain 
*(z + l) = (f 2 ^ 1 o F x )(z + l)+ v((* 2k ~ l o F A )(z + I)) (mod w^J 
= (* 2fc - 1 oF A )(z) + A(2a + 6) 



$(2 



(2«i + 1) 



+ A(2a + 6) 

+ ^ + 6 )(^TI) 



e- 1 (mod w^^J 



[(e x + e) + (e x - e) 

q 



+ i(2a + 6)L ± ^(2a + 6) - w Vl)W 

+ - ((2a + 6)L- bWt, 1)Vl ) (mod 
+ / (modw slil)1 ), 



(mod w Vl)t)1 ) 
(mod w„ liVl ) 



where we have used the fact that (e - * - + e) = (1,1), (e~ L — e) 
This completes the proof of theorem A. 



±(1, —1) and that p is odd. 

□ 



The set e of possible orbit codes is a subset of the product space 

2fe-l 

e e C{0,l,...,2^}x H {0,1,..., 2 Vj }. 

3=1 

Denoting by #L e the scaled co- volume of L e , namely 
(77) #L e = |AZ 2 /L e |, 

the total number of possible orbit codes is given by 



(78) 



|9 e | = #L e 



2^2 det (L, w Vl>Vl ] 



We note that although the lattice L e is independent of A (up to scaling), O e is not. 

In the next lemma, we identify the orbit codes which correspond to symmetric fixed points 
of Subsequently, in lemmas [13] and [bH we identify values of e for which the number of 
codes which satisfy the conditions of lemma [121 is independent of A. The proof of theorem B 
will then follow. 

Lemma 12. For any e E $ with vertex list (v±, . . . ,Vk), z G X e and sufficiently small A, the 
point z is a symmetric fixed point of $ if and only if its orbit code o~(z) = (<t_i, 0*1, . . . , o~2k—i) 
satisfies: 

(i) o-i = 01, (ii) 2a k = [\/eJ (mod 2\y/e\ + 1). 

Proof. It is a standard property of any reversible map F, with F = GoF- 1 oG, that it can 
be written as the composition of two involutions: 

G 2 = H 2 = Id. 



F = HoG, 



H = Fo G 



It follows that every symmetric periodic orbit of F intersects the union of the fixed spaces 
of these involutions, FixG U Fixff, at exactly two points (which coincide for period 1) [7J. 
Furthermore these two points are maximally separated in time, in the sense that if the minimal 
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period of the orbit is s, then the transit time from one to the other is approximately s/2. 
More precisely, if s is even, the orbit must intersect one of the following sets: 

(79) Fix G n F s / 2 (Fix G) or Fix if n F s / 2 (Fixif), 
whereas if s is odd, the orbit must intersect the set 

(80) FixG nF^ 1 )/ 2 (Fix if). 

For our map F\, we have already introduced the involution G and its fixed space, considered 
now as a subset of the rescaled lattice (AZ) 2 : 

G(x,y) = (y,x) FixG = {\(x,y) G (AZ) 2 : x = y}. 

A simple calculation shows that the involution if and its fixed space are given by: 

(81) H(x,y) = ([Xy\ - x,y) Fix if = {X(x,y) G (AZ) 2 : 2x = [Xy\}. 

Take e G $ and a point 2 G X e . Suppose that z is a symmetric fixed point of <£. As z is 
non-zero, the orbit of z intersects the set Fix G U Fix ff at exactly two points, and as z is a 
fixed point of these two points must occur within a single revolution. Hence we have: 

1 T (z) n (Fix G U Fix if )| = 2. 

We begin by considering which points in the return orbit O r {z) can lie in FixG. Since the 
domain X e lies in an 0(A)-neighbourhood of the positive half of the symmetry line FixG, we 
may have: 

z = $(z) G FixG. 

Equally, the orbit may intersect the negative half of the symmetry line, which occurs if: 

F 4 (z) = G(z) & F 2 (z) G FixG. 

Points in Fix H lie on disjoint vertical line segments of length one, in an 0(A)-neighbourhood 
of the y-axis. Recall that the polygon U(z) intersects the axes at vertices of type = \_\fe\ , 
and hence intersects the y-axis in the boxes: 

A),Lv^J and ^o-ciyij+i)- 
If L\/eJ is even, it follows that the relevant segment H\ of Fix H is given by: 

HI = {X(x,y) G (AZ) 2 : 2x = [Xy\ = Lv^J}, 

which lies in the positive half-plane. Similarly if [y/e\ is odd, the relevant segment Ht of 
FixH is given by: 

Ht = {X(x,y) G (AZ) 2 : 2x = [Xy\ = -([^\ + 1)}, 

which lies in the negative half-plane. 

Hence we see that the return orbit of z cannot intersect Fix H twice in a single revolution, 
and intersects FixG twice if and only if: 

z G FixG nr 2 (FixG). 

By (|79p . the latter implies that z is periodic with period 4, and we have already observed 
that there are no points with minimal period 4. Thus the only non-trivial possibility for a 
symmetric fixed point occurs when the return orbit of z intersects both FixG and Fix if. 
Conversely, equation (|80p ensures that this is also a sufficient condition. 
The proof now proceeds in two parts. 
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(i) O t (z) intersects FixG if and only if o--\ = oj. 

If z = X(x,y), then the property 0-\ = o\ is satisfied if and only if: 

x = y (mod 2v\ + 1) . 

If y = x then clearly z G FixG. The width of the strip X e , given by (|39p . ensures that the 
only other possibility is y = x + (2vi + 1), in which case: 

F 4 (z) = z + w„ 1)Ul 

= A(x + (2vi + l) J y-(2vi + l)) 
= A(y,x) = G(z). 

This corresponds to F 2 (z) G FixG. 

(wj O t (z) intersects FhtH if and only if2a k = LV^J (mod 2\^fe\ + 1). 

Instead of considering the sets and Ht. directly, we consider their images under G and 
F\, respectively, which lie in a neighbourhood of the x-axis: 

(82) G{H%) = {X(x,y) G (AZ) 2 : 2y = [Xx\ = Lv^J}, 

(83) F x (Ht) = {X(x,y) G (AZ) 2 : 2y = [-X(x + 1)J = -(Lv^J + 1)}- 

In l[83|) . we assume that A(|_\/eJ +l)/2 < 1, so that F\(w) = X(—l — y,x) for all w = X(x,y) G 
iff.. The orbit O t (z) intersects Fix// if and only if it intersects the relevant one of these sets, 
according to the parity of [v^J • 

The polygon H(z) intersects the x-axis at the kth vertex, where k is the length of the 
vertex list V(e). The return orbit O t (z) reaches the kth vertex at the point ^ k (z), given in 
the notation of (I47j) by 



^ k (z) = X 



+ x k ,Vk 



X 

where, by (i50l) . yf. = at is non-negative. Hence if [\/e\ is even, O t (z) intersects FixH if and 
only if: 

V k (z)eG(H e + ) a k =[^\/2. 
If \_y/e\ is odd, then O t (z) intersects FixH if and only if: 

F^ k (z)) G F\{Ht) & a k = -(Lv^J + l)/2 + (2 Wfc + 1) 

= (3Lv / ej+l)/2. 

The congruence 2cTfc = \^fe~\ (mod 2 [ v^J + 1) covers both of these cases, which completes the 
proof. D 



For all e G £ and sufficiently small A, the set X e — see equation (|39p — is non-empty and 
contains at least one element from every congruence class modulo L e . We now seek to identify 
the number of congruence classes whose orbit code satisfies the conditions of lemma [T2l 

As discussed in the proof of lemma [T2l the points z = X(x,y) G X e whose orbit code 
a(z) = (<t_i,<7i, . . . ,0"2fc-i) satisfies U-\ = o\ are precisely those satisfying 

y = x or y = x + {2v\ + l). 

All such points lie on one of two lines, parallel to the first generator L of the lattice L e . 
Furthermore all points on one line are congruent to those on the other, as they are connected 
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by the second generator (L — w VlVl )/2. Hence the number of points satisfying this condition 
modulo L e is 

(84) |L| = = _£_, 

where we have used ([75]) . 

It remains to determine what fraction of the orbit codes with o~-\ = o~\ satisfy the second 
condition of lemma [12] We do this by identifying values of e for which all possible values of 
cr fc occur with equal frequency, independently of A. 

Clearly if k = 1, i.e., if a polygon class has just one vertex in the first octant (e = — a 
square), then the points z = X(x,y) G X e with ay. = o\ = a* satisfy 

x = y = a* (mod 2vi + 1) , 

for any given a* G {0, 1, . . . , 2v\}. Such points form a fraction 



{2v x + l) 2 

of all points modulo L e . Hence all possible values of occur with equal frequency modulo 
L e . More generally if = v\, i.e., if all vertices of the polygon class have the same type 
(e = 0,2,8), then the same applies. This follows from the fact that, for any congruence class 
of L e , the map Oj t-t (Jj+i is a permutation of the set {0, 1, ... , 2vj} whenever Vj = Vj+i, as 
we saw in case 1 of the proof of lemma [10] 

The following lemma deals with the case that a polygon class has two or more distinct 
vertex types. 

Lemma 13. Let e G $ . Suppose that the vertex list (t>i,t>2, • • • ,Vk) of the associated polygon 
class has at least two distinct entries and satisfies 

(85) gcd(2^ (0 + l,p ty _ 1) ) = l, 

where (i(i)) i=1 is the sequence of distinct vertex types defined in (T53]). Then for every z G X e , 
all 1 < j < l(1), all a* £ {0, 1, . . . , 2vt}, and all sufficiently small X, the number of points in 
the set (z + L|)/L e whose orbit code has kth entry a* is 

(86) 2^TT |L 5/ Le| ' 



Proof. Suppose that the coprimality condition (|85j) holds. Let z G X e have orbit code 
a(z) = (<7_i, <7i, . . . , 0"2fc-i) and let the pair (xj,yj) be defined as in equation ([17]) . 

The statement (1861) is equivalent to all possible values of cifc occurring with equal frequency 
among points in z + modulo L e . Hence it suffices to prove that (|86j) holds for j = l(1) — 1, 
since all cylinder sets of L| with index i < j can be written as a union of cylinder sets of L^. 

We let j = l(1) — 1, so that 

v 3 + l = v l{1) = v k, 

and consider the congruence class of z modulo L^. Since z + is a cylinder set in the sense 
of lemma [TU1 the orbit codes of all points z G z + match up to the jth entry o~j . Let the 
(j + l)th entry of the orbit code of such a z be We wish to show that all possible values 

of crj + i occur with equal frequency. 

By construction Vj / Vj+i, so the possible values of are determined by case 2 of the 
proof of lemma HU] In the course of the proof, we saw that the occurrence of points z with 
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some fixed value of CTj+i correspond to solutions of an integer equation, given in the case 
where y is the non-integer coordinate of the jth vertex by equation (|68p . (A similar equation 
holds when x is the non-integer coordinate.) Each solution (2a + b, t) G Z x Kl determines the 
module coordinates (a, 6) of z — z in L| and the transit time t of 5 from the jth vertex to the 
(j + l)th. 

Solutions of ([68]) occur for all values of <5j+i satisfying the congruence ([69]) . and the condi- 
tion that A be sufficiently small ensures that at least one such solution is realised by a point 
z G X e . By construction, each distinct value of which has a solution defines a unique 
point in z + IL| modulo L| +1 , which is isomorphic to the module DJ-/Lj +1 . However due to 
the coprimality condition (|85|) . the modulus of the congruence (j69|) is unity. Hence points z 
occur for all possible values of (Xj+i, and each is unique in z + QJ? modulo D-^+i- Furthermore, 
by ([62]) . the lattices L^ +1 and L e are equal, hence all possible values of occur with equal 
frequency in z + modulo L e . 

If j + 1 = = k then this completes the proof. If l(1) < k, take i in the range t{l) < i < k. 
By the definition of as the index of the last distinct vertex type, we have Vi = Vi + \ = Vk- 
As discussed above, the map <7j h-> a,i + \ is a permutation of the set {0, 1, ... , 2vi} whenever 
Vi = Vi + \. Hence the equal frequency of the possible values of cr, implies that of and the 
result follows. D 

In the previous section (equation (jlHj) ). we defined the jth entry cij of the orbit code via the 
congruence aj = yj (mod 2vj + 1), where y is the integer coordinate of the relevant vertex, 
and the pair (xj,yj) is defined by equation (|T7j) . Similarly, we define the sequence 7(2) such 
that its jth entry jj satisfies: 



(87) 7j = Xj I mod 



2vj + 1 



where x is the non-integer coordinate of the vertex. It follows from corollary 1111 that . for any 
1 < j < 2k — 1, and any two points z, z G X e : 

z = z (mod L e ) <^> (<Tj, 7,) = (^j,7j)- 

We use 7(21) to show that, among points with a-\ = u\ and for each j in 1 < j < 2k — 1, all 
possible values of <7j G {0, 1, . . . ,2t>j} occur with equal frequency modulo L e , independently 
of A. 

Lemma 14. Let e G and suppose that the vertex list (vi, V2, ■ ■ ■ , ffc) 0/ i/ie associated 
polygon class is such that 2vi + 1 is coprime to 2vj + 1 for all other vertex types Vj : 

(88) gcd(2ui + 1, 2vj + 1) = 1 2 < j < k, vj / v x . 

Then for sufficiently small A, for all j in 1 < j < 2k — 1, and all a* G {0, 1, . . . ,2vj}, the 
number nj of points z G X e modulo L e whose orbit code (cr_i, a\, . . . , 0"2jfc_i) /ias <r_i = <ti 
and aj = a* is given by: 



(89) 



(2wi + l) (2uj + l)' 



Proof. We use induction on j. Consider points 2 whose orbit code satisfies cr_i = o"i and 
has jth value <Tj, for some arbitrary aj G {0, 1, ... , 2vj) and j G {1, . . . , 2k — 1}. Let the 
sequence 7(2) be denoted (7-1,71, . . . ,72fc-i)- Our induction hypotheses are that: 
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(i) equation (|89[) holds, coprimality condition (|88p ensuring that rij is a natural number; 

(ii) for each residue r £ {0, 1, . . . , rij — 1} modulo rij, there is a unique z modulo L e satisfying 

7j = r (mod n^) . 

The base case is j = 1. The points z = A(x, y) E X e with cr_i = o"i for some fixed value of 
o\ satisfy: 

x = y = a\ (mod 2v\ + 1) . 

Such points are congruent modulo (\(2v\ + 1)Z) 2 , hence the number of such points modulo 
L e is 

#L e g 
(2 Vl + l) 2 (2 Wl + l) 2 ni " 
By lemma [TOl if z is one such point, then any other point z reaches the first vertex at 

= V(z) + Xspie 

for some s 6 Z, where pi = 2ui + 1 and e is the unit vector in the non- integer coordinate 
direction. Furthermore, if 7(2) = (7-1,71, . . . ,72fc-i), the value of 71 is related to 71 by 

7i = 7i + s(2v\ + 1) (mod (2v\ + 1) n{) . 

By corollary [TTJ z and z are congruent modulo L e if and only if 71 = 71 . Thus the n\ distinct 
points modulo L e correspond to distinct values of s modulo n\. Furthermore, if we consider 
the value of 71 modulo m, we have: 

7i = 7i + s{2v\ + 1) (mod ni) . 

Now 2t>i + 1 is coprime to the modulus, as by the coprimality condition (188p and the con- 
struction (155p of g, two is the highest power of {2v\ + 1) that divides q. It follows that each 
distinct value of 71 is distinct modulo m. This completes the base case. 

To proceed with the inductive step, we suppose that the above hypotheses hold for some 
j G {1, . . . , k — 1}. In the proof of lemma [TOl we used equation (I52p to describe the behaviour 
of points as they move from one vertex to the next in two cases. The first case occurs when 
Vj = Uj+i, so that the jth and (j + l)th vertices lie on parallel lines, and rij = rij + \. In this 
case, the value of crj+i is determined uniquely by the value of cr,-. In particular, we saw that 
if the jth. vertex lies on y = n and the (j + l)th vertex lies on y = n — 1, then <7j+i and <7j 
are related by equation (p7|) . 

We can use the same methods, considering this time the non-integer component of equation 
([52]) . to show that 7 J+ i is determined by the pair (<Tj,7j) via: 

7 J+ i = 7j + €j + t(2n - 1) (mod (2v\ + 1) rij) , 

where ej = €j(aj), and t = t(aj) is the transit time between vertices. The one-to-one rela- 
tionship between aj and (Jj+i, ensures that there are %-+i = rij points with cr_i = <ri that 
achieve any given value of &j+i at the (j + l)th vertex. Then for any given value of aj, 
which uniquely determines (Jj+i, the above congruence establishes a one-to-one relationship 
between jj and Jj+i modulo (2vi + l)rij. Because this bijection is a translation, it also 
holds modulo rij. In other words, there is a skew-product map of residue classes modulo rij: 
(<7j,7j) 1 y (oj+i, 7j'+i))- This completes the inductive step for the first case. 

In the second case, where Vj 7^ Vj+i, the jth and (j + l)th vertices lie on perpendicular 
lines. Again referring to the proof of lemma [TOl taking equation (|68|) modulo 2vj + \ + 1 gives 
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the following expression for Oj+i in terms of the pair (ajjjj): 



A 



+ (Tj+l 



+ 7j + ej (mod 2vj + i + 1) 



Here we were able to replace Xj with jj as, by the construction (|55p of q, 2vj+% + 1 is a divisor 
of the modulus q/{2vj + 1) = {2v\ + l)n,- which defines jj. If the coprimality condition 
([88]) holds, then 2vj + i + 1 also divides rij. Hence for any given pair (aj,aj+%), there are 
rij/(2vjJ r i + 1) values of 7j modulo nj for which the following congruence is satisfied: 

- Vj + 1 



(90) 



7? 



A 



A 



6j + s(2vj + i + 1) (mod no- 



where s € Z. The total number of points with any given value of Cj+i is thus: 



(2vj + 1) x 



■j+i, 



2v j+1 + 1 

which completes the inductive step for hypothesis (i). 

Taking the second component of equation ([52]) modulo rij+i gives an expression for Jj+i 
in terms of the pair (<jj , jj ) : 

71—1 

+ <Tj + i(2«j + 1) - 



(91) 



7j+l 



A 



(mod rij+i) 



where t = t(aj,jj). For a given pair {<jj, <Jj+i), t is given by equation ([S5]) . Hence taking 
equation ([68]) modulo rij, a multiple of 2vj + i + 1, and using the expression ([90]) for jj, it 
follows that the values of t satisfy 

K^ + 1)/A1 +a j+1 -\ Vj /X\ 



t + l = 



7j 



2v 



+ 1 



mod 



7 J ; 



2?> 



'i+i 



+ 1 



= -s (mod rij/ (2vj + \ + 1)) , 

where s G Z. Thus £ takes all values modulo rij/(2vj + \ + 1) = rij+i/(2vj + 1). In turn, 
equation ([9~T]) implies that 'Jj+i takes all values satisfying: 

n — 1 



7j+i = 



A 



(mod 2vj + 1) 



Applying this argument for all values Oj £ {0,1,... ,2fj}, we get a complete residue class 
modulo nj + i as required. This completes the inductive step for hypothesis (ii) and the result 
follows from hypothesis (i) for j = k. d 

Proof of theorem IB] 

Let e £ if be given, let (ui, . . . ,ffc) be the corresponding vertex list, and let <r* be the 
unique element of the set {0, 1, ... , 2vk} that satisfies 

2a*=[^\ (mod2Lv^J +!)• 

For e = 0, 2, 8 all elements of the vertex list are the same. This case is dealt with by the 
discussion preceding lemma [13l Thus we assume that the vertex list contains at least two 
distinct elements. 

Suppose first that 2v\ + 1 is coprime to 2vj + 1 for all Vj ^ v±. Then for sufficiently small A, 
lemma [HI states that the number of points in X e n FixG modulo L e whose orbit code has 
kth. entry a* is given by equation ([89]) . with j = k. Furthermore, since all points z £ X e whose 
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orbit code (<r_i, 01, . . . , (T2k—l) satisfies cr_i = o\ are congruent to some point in X e n Fix G, 
it follows that n k is the number of points in X e modulo L e satisfying the conditions (i) and 
(ii) of lemma [12j Therefore the number of symmetric fixed points of $ in X e modulo L e is 
independent of A and given by n k - 

Similarly if 2v k + 1 is coprime to 2vj + 1 for all Vj ^ v k , then 2v k + 1 is coprime to q L n-2), 
given in closed form by (|6ip . If follows that the condition (|85p of lemma [T3l holds, since the 
recursive expression ([57|) for g t (/_i) gives us that: 

Qi(Z-l) 

Pi(J-l) " 



_ lcm((2z; t (;_ 1) + 1)(2^_ 2) + l),g t (;_ 2 )) 
2u t( ,_ 1) + 1 

_ (2^-2) + 1) 

gcd((2u i(z _ 1) + l)(2u t(z _ 2 ) + l),&(7-2))' 

Applying lemma [Till for j = 1, we have that for every cylinder set of Lf , the number of points 
modulo L e in the cylinder set whose orbit code has kth entry a* is given by 

1 |L e /L e, 1 #^ 



2v k + V u 1 2v k + 1 (2v! + l) 2 2«i + l" 

There are 2v± + 1 cylinder sets of Lf whose associated orbit code satisfies cr_i = o\. Hence, as 
before, the number of points in X e modulo L e satisfying the conditions (i) and (ii) of lemma 
E2]is 

(2vt + 1) x fc = n k , 
2v\ + 1 

and the number of symmetric fixed points of $ in X e modulo L e follows. This completes the 
proof of the first statement. 

We have shown that for sufficiently small A, and if (142ft holds, then the fraction of symmetric 
fixed points of in each fundamental domain of L e is 



(2«i + l)(2« fc + l) {2[^j2\ +l)(2Lv^J +1)' 

where we have used equations (|22p and (|23p . respectively. It remains to show that the density 
8(e, A) of symmetric fixed points in X e converges to this fraction as A — > 0. 

By equation ([39]) . the domain X e is a subset of the lattice (AZ) 2 , bounded by a rectangle 
lying parallel to the symmetry line FixG. Similarly, a fundamental domain of the lattice L e 
is constructed by taking the points bounded by the parallelogram £2, given by 

n = {aL + ^(L - w VuVl ) : a,/?e[0,l)}, 

where the generator L is also parallel to the symmetry line. These parallelograms tile the 
plane by translation: £1 + L e = R 2 . 

The width of X e (taken in the direction perpendicular to FixG) is exactly twice that of £1, 
independently of A, as shown in Figure [9l The number of parallelograms which fit lengthwise 
into X e , however, goes to infinity as A goes to zero. In particular, the number of parallelogram 
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which can be contained in the interior of the rectangle bounding X e is given by 

\I e \ 



A|L| 

where |J-f e |/A|L|J is the number of times that the vector L fits lengthways into the rectangle, 
and we subtract 8 for the parallelograms which intersect the boundary. Each of these par- 
allelograms contains a complete fundamental domain of L e , and their contribution to <5(e, A) 
dominates in the limit A — > 0. 
Explicitly, we have: 



5(e,\) 



#L e / l |_| J C |/A|L| 

JT* [ {2 Vl + l)(2v k + 1) 



+ 0(1) 



{2v l + l)\~F\/\ q 



2(2^ + 1) |/ e |/A 



(2«i + l){2v k + 1) 



(2«i + l)(2« fc + 1) 
+ 0(A), 



+ 0(1) 



as A — > 0. 



□ 
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